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1 Introduction 



Physics of domain walls and vortices has been a rather interesting field of 
experimental as well as theoretical research for quite a long time ^ . Example 
of a recent hot topic is production of such soft solitonic objects in rapid 
phase transitions — testing and refining a theoretical description proposed 
by Kibble and Zurek [0]. Theoretical analysis of dynamics of domain 
walls and vortices is relatively difficult because pertinent field equations are 
nonlinear, and the most interesting solutions do not belong to weak field 
sector. 

Among problems which have been discussed in literature is time evolution 
of single curved domain wall or vortex. It is commonly regarded as accessible 
only by a numerical analysis. Actually, there exist also analytical approaches 
which yield (an approximate) description of the time evolution: the classical 
effective action (CEA) method which has been developed in a series of pa- 
pers starting from [0,0, examples of more recent works are 0, and 
a version of Hilbert — Chapman — Enskog method which we call the im- 
proved expansion in width (JEW). This latter approach has been developed 
in papers ||^, with an inspiration coming from The two methods have 



been outlined and compared in paper |TT|. These two analytical approaches 



are not simple, but neither is the purely numerical approach — this is just a 
reflection of the fact that dynamics of the curved domain wall or a vortex is 
nontrivial due to nonlinearity, many spatial dimensions, and many modes of 
involved flelds. The numerical and analytical approaches should be regarded 
as equally important and complementary sources of information about the 
dynamics. 

In the present paper we apply lEW method to a curved domain wall in a 
nonrelativistic dissipative system. Time evolution is governed by a diffusion 
type equation for which no simple action functional exists. Therefore, it 
is not clear how CEA scheme could be applied in this case, while, as it 
turns out, lEW method works quite well. Because our second goal is a 
presentation of the method, we consider a relatively simple system with scalar 
order parameter. An application to domain walls in nematic liquid crystals 
we will present elsewhere |T^ . lEW method has also been applied to a vortex 
line, see |T^. Work on application to a disclination line in a nematic liquid 
crystal is in progress [|14| . 



In our opinion lEW method has several attractive features, e.g., it com- 
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bines the old and elegant subject of differential geometry of surfaces in 3- 
dimensional space with nonlinear dynamics of the curved domain wall. An- 
other interesting aspect is that lEW scheme relates properties of the domain 
wall to properties of one-dimensional kink. In this sense JEW method embod- 
ies the idea that the curved domain wall can be regarded as three-dimensional 
embedding of the one-dimensional kink. 

The expansion in width is based on the idea that transverse profile of 
the curved domain wall considered in suitable coordinates (which are called 
comoving coordinates) differs from transverse profile of a planar domain wall 
by small corrections which are due to curvature of the domain wall, and that 
these corrections can be calculated perturbatively. There is a condition for 
applicability of such a perturbative scheme: the two main curvature radia of 
the domain wall should be much larger than its width. As we shall see below, 
turning that idea into a concrete calculational scheme requires some work, 
but that should be expected in any approach which tackles generic curved 
domain walls. 

The plan of our paper is as follows. In Section 2 we introduce the co- 
moving coordinate system. Section 3 is devoted to the presentation of the 
perturbative expansion. Several remarks are collected in Section 4. 



2 The comoving coordinates 

We shall seek the curved domain wall solutions of the following equation 
where the free energy F has the form 

Here $ is a real scalar order parameter; A, v and 7 are positive constants; 
and {x°'), a = 1,2, 3, are Cartesian coordinates in the usual space. From 
(1) and (2) we obtain the following equation for the rescaled dimensionless 
order parameter (p = ^/v 

7^ = A</.-2A^;W-l), (3) 
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where A = dada- The domain wall solutions of Eq.(3) smoothly interpolate 
between = —1 on one side of the wall and = +1 on the other side. 
Canonical example of such a solution is given by the formula 

3 

0o(a^^) =tanh ^ ^ , (4) 

where Iq"^ = 4Aw^ and a is an arbitrary constant. The presence of the constant 
a is due to translational invariance of Eq.(3). This particular solution rep- 
resents planar static domain wall located at the plane = a. Such domain 
wall is homogeneous along that plane. Its transverse profile is parametrized 
by x^. Width of the wall is approximately equal to lo, in the sense that for 
— a| ^ /q values of differ from +1 or -1 by exponentially small terms. 
The form (4) of (poix^) coincides with one-dimensional static kink present in 
one-dimensional version of the model defined by (1) and (2). The domain 
wall can be regarded as embedding of that kink in the 3-dimensional space 

Notice that somewhere inside any domain wall there is a surface on which 
vanishes. For example, in the case of planar domain wall (4) 0o = for 
x^ = a. Such surface is called the core of the domain wall. 

The first step in our construction of the perturbative scheme consists in 
introducing special coordinates comoving with the domain wall. One coor- 
dinate, let say ^, parametrizes direction perpendicular to the domain wall, 
two other coordinates (cr^, cr^) parametrize the domain wall regarded as a 
surface in the space. The comoving coordinates have been proposed in 
the context of CEA method, We introduce one important modification: 
an auxiliary surface S, which is present in the definition of the comoving 
coordinates, is apriori independent of the domain wall. In literature on CEA 
method it is defined directly in terms of the domain wall — the most popular 
choice is that S coincides with the core. It turns out that the latter choice in 
general is not compatible with certain consistency conditions which appear 
in our approximation scheme. Transformations to comoving coordinates in 
the cases of relativistic domain walls and vortices in Minkowski space-time 
can be found in 0, |T^, respectively. Below we introduce such coordinates 



for the nonrelativistic domain wall moving in the R space. 

After all these remarks let us finally define the comoving coordinates. We 
consider a smooth, closed or infinite surface S in the usual R^ space. It is 
close to the core and its shape roughly gives the shape of the domain wall. 
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In particular one may assume that S coincides with the core at certain time 
to- Points of S are given by X{a^,t), where a* {i = 1,2) are two intrinsic 
coordinates on S, and t denotes the time — we allow for motion of S in the 
space. The vectors X^, k = 1,2, are tangent to S at the point X{a^,t) []. 
They are linearly independent, but not necessarily orthogonal to each other. 
At each point X{a\t) of S we also introduce a third vector 'p{a\t) which is 
perpendicular to S, that is 

pXfc = 0. 

We assume that p has unit length, jp' = 1. The three vectors {X^k,p) form a 
local basis at the point X(cr*,t) of 5*. With this basis given at each point of 
5*, we introduce geometric characteristics of S: induced metric tensor on S 

Qik = X^iX^ki 

and the extrinsic curvature coefficients 

Kii = pX^u 

{i, k,l = 1,2). They appear in Gauss- Weingarten formulas 

X,ik = K,kP + r[,X;, = -fKikX^i. (5) 

Here the matrix {g^^) is by definition the inverse of the matrix {gki)-, i-e. 
g^^Qki = ^l) and T\f, are Christoffel symbols constructed from the metric Qik- 
The two by two matrix {Kik) is symmetric. Two eigenvalues ki,k2 of the 
matrix (A'j), where i^j = g'^^Kij, are called extrinsic curvatures of 5* at the 

point X. The two main curvature radia are defined as Ri = 1/ki. In general 
they vary along S and with time. 

The comoving coordinates {a^,a^,^) at the time t are introduced by the 
following formula 

x = Xia\t)+^pia\t). (6) 

Here x = = 1,2,3, are the usual Cartesian coordinates in the space 

R^. ^ is the coordinate in the direction perpendicular to S. Notice that this 
direction has very simple parametrization — the r.h.s. of formula (6) is a 
linear function of ^. We will use a compact notation for the comoving coor- 
dinates: ((T^,(T^,^) = (ex"), with a=l, 2, 3 and = ^. The coordinates (a") 

^We use the compact notation f^k = df/da'^. 
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are just a special case of curvilinear coordinates in R^. The corresponding 
metric tensor Gap in has the following components: 



G33 — 1, Gsk — Gk3 — 0, Gik — N^girN^, 



where 




where G = detiGap), g — det{gap), and N — det{Nl) is given by the 
following formula 



We see that dependence of Gap on the transverse coordinate ^ is explicit, 
and that a^,a^ appear through the tensors g^k, which characterise the 
surface S. 

In general the coordinates (cr") have certain finite region of validity. In 

particular, the range of ^ is given by the smallest positive ^o(c\ t) for which 
= 0. It is clear that such ,^0 increases with decreasing extrinsic curvature 
coefficients K^, reaching infinity for the planar domain wall. We assume that 
the surface S (and the domain wall) is smooth enough, so that outside of 
that region there are only exponentially small tails of the domain wall which 
give negligible contributions to physical characteristics of the domain wall. 

The comoving coordinates are utilised to write Eq.(3) in a form suitable 
for calculating the curvature corrections. Laplacian in the new coordi- 
nates has the form 
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Components G"^ of the inverse metric tensor have the form 



^ ^ Q3k ^ Qk3 ^ Qik ^ ^^-l^ji^rt^^-l^fe^ 



where 



{N-y, = - {{1 - ^Kl)6i. + ^K) . 




6 



The time derivative on the l.h.s. of Eq.(3) is taken under the condition that 
all a;° are constant. It is convenient to use time derivative taken at constant 
cr"^. The two derivatives are related by the formula 

d . d . da^ d 

where 

the dots stand for d/dt\^i. The final step consists in rescaling the transverse 
variable ^ 

^ = 2/os. 

The dimensionless variable s measures the distance from the surface S in 
the unit 21q. Equation (3) transformed to the comoving coordinates with ^ 
rescaled as above has the following form 

which is convenient for construction of the expansion in width. 



3 The expansion in the width 

We seek domain wall solutions of Eq.(7) in the form of expansion with respect 
to Zo, that is 

= 00 + + ^0</'2 + ■- (8) 

Inserting formula (8) in Eq.(7) and keeping only terms of the lowest order 
(~ Iq) we obtain the following equation 

i^ + ^o-^^o. (9) 

It has the well-known kink solutions 

00 = tanh(s — sq), 
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which formally have the same form as the planar domain walls (4). In the 
remaining part of the paper we shall calculate curvature corrections to the 
simplest solution 

0o = tanhs. (10) 

Notice that 0o interpolates between the vacuum solutions ±1. Therefore, the 
corrections k > 1, should vanish in the limits s ±oo. 

Equations for the corrections (f)k, k > 1, are obtained by expanding the 
both sides of Eq.(7) and equating terms proportional to Iq. They can be 
written in the form 

Uk = fk: (11) 

where 

and fk depends on the lower order contributions 4>i, I < k. Straightforward 
calculations give 

h = d,MK-iP^), (13) 

/2 = 3000? + 2sds(t>oK)Ki + dsMK - IPX), (14) 



h = 27(9*01 - g^'X^rXdM + 6000102 + 0? 
+2sds(t>iKiKi - 2s29,0oX; [{Klf - SKjKf) 

-j^d.i^g^'^dM + dsMK - IP^)^ (15) 



and 



U = 27(5*02 - 2sg''^'pX^A<pi) - 27(?^-^X fcX(9,02 + 2sKp,<P,) 
+30002 + 6000103 + 30202 + 2sd,(t)2KiKi 
-As^ds4>o {{K;Y - {KlK^.f - 2{K;fK]Kl) 
-2s^ds(t>,Kl {{KlY - SKjKi) - j^dji^g^'^dM 
-f^d,{^Ki^dk^i) + Asg^%K;dk4>i + dsHK - IPX), (16) 
where dt — d/dt, di — d/da^. 
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Notice that all Eqs.(ll) for 0^ are linear. The only nonlinear equation in 
our perturbative scheme is the zeroth order equation (9). 

Now comes the crucial point: operator L has a zero-mode, that is a 
normalizable function iI)q{s) such that 

L^o = 0. 

This function can be obtained by differentiating (t)Q{x^) given by formula (4) 
with respect to a/2/o and putting a = 0, 

% = (17) 

cosh s 

This zero- mode owes its existence to invariance of Eq.(3) with respect to 
spatial translations, therefore it is often called the translational zero-mode. 
Let us multiply Eqs.(ll) by ipo{s) and integrate over s. Because 



we obtain the consistency (or integrability) conditions 

/oo 
dsMs)fkis) = 0. (18) 
'OO 

The operator L appears also in the expansion in width for relativistic do- 
main walls 0. Using standard methods |1^,0 one can obtain the following 
formula for vanishing in the limits s —>■ ±oo solutions (pk of Eqs.(ll): 

<Pk = G[fk] + Ck{a\T)Ms), (19) 

where 

G[fk] = -2Ms) r dxMx)fk{x) + 2Ms) r dxMx)fkix). (20) 



Here ipo{s) is the zero-mode (17) and 



V'i(s) = isinh(2s) + ^tanhs + ^^^ (21) 
o 8 8 cosh s 



is the other solution of the homogeneous equation 

Lip = 0. 



9 



The second term on the r.h.s. of formula (19) obeys the homogeneous equa- 
tion = 0. It vanishes when s ±00. 

One can worry that (pk, k > 1, given by formulas (19), (20) do not vanish 
when s ±00 because the second term on the r.h.s. of formula (20) is 
proportional to ipi which exponentially increases in the limits s ±00. 
However, the integrals 

dxipofk 



vanish in that limit, see the consistency conditions (18). Moreover, qualita- 
tive analysis of Eq.(7) shows that fk ~ (polynomial in s) x exp(— 2|s|) for 
large |s|, hence those integrals behave like (polynomial in s) x exp(— 4|s|) for 
large \s\ ensuring that (pk exponentially vanish when |s| — > 00. 

We have explicitly solved Eqs.(ll). The solutions (19) contain as yet 
arbitrary functions Ck{(j\t), and also X{a\t) giving points of the comoving 
surface S — K^, go. follow from X. It turns out that the conditions (18) 
are so restrictive that they essentially fix these functions. The consistency 
condition with k = 1 is equivalent to 

jpJl = K;, (22) 

where we have used formulas (10), (13) and (17). Thus we have obtained 
equation for X. It is of the same type as AUen-Cahn equation |T6[, but in our 



approach it describes motion of the auxiliary surface S only. Equation (22) 
should be compared with Nambu-Goto equation for a relativistic membrane 
obtained in the relativistically invariant version of our model . We shall see 
below that the remaining consistency conditions do not give more restrictions 
for X at least up to the fourth order — they can be saturated by the functions 
Ck{cr\ t). We expect that this is true to all orders but we have not attempted 
to provide a proof. 

Let us now proceed with the discussion of the perturbative corrections: 
Taking into account the condition (22) we have /i = 0. Therefore 

^ (23) 

cosh s 

Equation (11) with k = 1 does not provide more information. 



10 



The second order contribution 02 is given by formula (19) with k = 2. 
Using the results (22, 23) from the first order we obtain the following formula 

02 = Ms)C!{a\t)+Us)KiKi + (24) 

cosh s 

where 

sinhs 
cosh^ s ' 

, / N -4 /•« xiJi(x) , , , . r , ^ 

V'3 s) = — -2- / dx — -2— + 4V^i s / dx — -4—. 
cosh s Jo cosh X J-00 cosh x 

■03 (s) can easily be evaluated, e. g., numerically. Also the higher order 
corrections involve only rather simple integrals of elementary functions. 

The consistency condition (18) with k = 2 does not give any restrictions 
— it can be reduced to the identity = 0. More interesting is the next 
condition, that is the one with k = 3. It can be written in the form of the 
following inhomogeneous equation for Ci{a\ t) 

7(9tCi - g'''X,r^^kC^) - j=dj{^g^'^dkC,) - K^KiC 

= - 1)^; {(Kir - 3K]Kl) . (25) 

This equation determines Ci provided that we fix initial data for it. The 
consistency condition coming from the fourth order (/c = 4) is equivalent to 
the following homogeneous equation for C2 

l{dtC2 - g'''-X,JdkC2) - -^d,{^g^%C2) - K}KiC2 = 0. (26) 

The perturbative scheme presented above is not quite straightforward. 
Therefore we would like to add several explanations. The formulas presented 
above give a whole family of domain wall solutions. To obtain one concrete 
domain wall solution we have to choose initial position of the auxiliary surface 
S. Its positions at later times are determined from Eq.(22). We also have 
to fix initial values of the functions Ci,C2 and to find the corresponding 
solutions of Eqs. (25), (26). The approximate domain wall solution (p is 
given by formulas (8), (19), (20). Notice that we are not allowed to choose 
the initial profile of the domain wall arbitrarily because the dependence on 



V'2(s) 
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the transverse coordinate s is explicitly given by these formulas. Any choice 
of the initial data gives an approximate domain wall solution. Of course 
such a choice should not lead to large perturbative corrections at least in 
certain finite time interval. Therefore one should require that at the initial 
time IqCi ^ 1,IqC2 ^ ^ 1- The domain wall is located close to the 

surface S because for large |s| the perturbative contributions vanish and the 
leading term tanhs is close to one of the vacuum values ±1. 

The presented formalism is invariant with respect to changes of coordi- 
nates a^. (7^ on S. In particular, in a vicinity of any point X of S we can 
choose the coordinates in such a way that gik — Sik at X. In these coordinates 
Eq.(22) has the form 

7^ = -^ + -^, (27) 
rCi 1X2 

where v is the velocity in the direction p perpendicular to S and R2 are 
the main curvature radia of S at the point X. 

Let us present a simple example: take S to be a sphere of radius R. Then 
R^ = R^ = -R{t), t; = i? and Eq.(27) gives 



R{t) = ^Rl-^{t-to), 

where i?o is the initial radius. Our approximate formulas are expected to be 
meaningful as long as R{t)/lo » 1. Equations (25), (26) reduce to 

In the last equation 6', are the usual spherical coordinates on S (we apologize 
for using letter also in this meaning) . If we take the simplest initial data at 
t — to, namely Ci — C2 — 0, then the last equation imphes that C2 — also 
for t > to while Ci > 0. The Cartesian coordinate frame is located at the 
center of the sphere S and p is the outward normal to 5"; s — {r — R{t))/2lo, 
where r is the radial coordinate in R^. 
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4 Remarks 



We would like to make several general remarks about the expansion in width. 

1. We have used Zq as a formal expansion parameter. It is a dimensionful 
quantity, hence it is hard to say whether its value is small or large. What 
really matters is smallness of the corrections lo(pi, Iq4'2- This is the case if 
IqCi < 1, IIC2 < 1 and IqK] < 1, as seen from formulas (8), (19) and (20). 

2. Notice that an assumption that S coincides with the core for all times in 
general would not be compatible with the expansion in width. If we assume 
that Ci = = C2 at certain initial time toi Eq.(25) implies that Ci 7^ at 
later times (unless the r.h.s. of it happens to vanish). Then, it follows from 
formulas (8), (10) and (23) that does not vanish at s = 0, that is on S. 

3. The question of convergence of the expansion (8) has not been analysed. 
Actually we think that the expansion can turn out to be convergent, in spite 
of the fact that more frequent in field theory are asymptotic expansions. 
Moreover, this problem seems to be within the reach of the present day 
mathematical techniques. 

4. Finally, we would like to stress that we have abandoned effects which come 
from perturbations of the exponential tails of the domain wall. For example, 
if we have a domain wall in the form of cylinder with very large radius and 
small height (and with rounded edges), then the top and bottom parts are 
flat, and according to Eq.(27) they do not move. In our approximation the 
cylinder shrinks from the sides where the mean curvature I/R1 + I/R2 does 
not vanish. Now, in reality the top and bottom parts interact with each 
other. This interaction is very small only if the two flat parts are far from 
each other. We have neglected it altogether assuming the tanh s asymptotics 
for large s. Thus, our approximate solution takes into account only the effects 
of curvature. 
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